Universality of Sudden Death of Entanglement 
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We present a constructive argument to demonstrate the universality of the sudden death of entan- 
glement in the case of two non-interacting qubits, each of which generically coupled to independent 
Markovian environments at zero temperature. Conditions for the occurrence of the abrupt dis- 
appearance of entanglement are determined and, most importantly, shown to be almost always 
satisfied: Dynamical models for which the sudden death of entanglement does not occur are seen to 
form a highly idealized zero-measure subset within the set of all possible quantum dynamics. 
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I. INTRODUCTION 

The entanglement of quantum states is one of the 
most fundamental aspects of quantum mechanics and 
is crucial for countless processes of quantum informa- 
tion and computation ||, ^, || ||, ^. Recently, the 
study of entangled states perturbed by local environ- 
mental noise has become a subject of increasing interest. 
Unlike in decoherence processes 0, there are some en- 
tangled states which, depending on their couplings to the 
environment, become separable within a finite time inter- 
val § |, |l0l m, tl|, 111, 111- Such states manifest 
what is called the sudden death of entanglement (SDE). 

The incidence of SDE has been studied for a large class 
of initial entangled qubit states, and in the separate cases 
of amplitude-damping, depolarizing, dephasing, bit flip- 
ping, and bit-phase flipping |, |, jg, |l| ^. 

However, of all possible coupling agents inducing the 
noisy dynamics, the ones associated with these sepa- 
rate cases constitute very particular choices and, in many 
cases, not realistic ones. 

As an example, consider the case where a photon, in 
a linear combination of zero- and one-photon states, tra- 
verses an optical medium and is scattered by an atom. In 
such a process, as it is well known, dispersion is positively 
correlated with absorption, implying that we cannot re- 
alistically describe the environmental noise as amplitude- 
damping alone; we must also consider dephasing. 

Motivated by practical constraints such as above, and 
also by the theoretically large set of available coupling 
agents between the system and the environment, we are 
led to ask whether SDE is a general feature of quantum 
dynamics, or else a peculiarity of certain models. 

In this paper, this problem is approached by consid- 
ering generalized error agents in the case of two non- 
interacting qubits, each of which coupled to its own in- 
dependent Markovian environment at zero temperature. 
On top of deriving conditions for the occurrence of SDE, 
we show that the set of coupling agents for which SDE 
may not occur has zero measure with respect to the set of 



all possible couplings, thus proving that — in the present 
circumstances — SDE is a universal phenomenon. 

The paper is organized as follows. In Sec. |l| we intro- 
duce the qubit interaction Hamiltonian for general error 
agents, write and solve the Born-Markov master equa- 
tion for the dynamics of a single qubit, and review a 
procedure to construct the dynamics of two independent 
qubits from the single qubit dynamics. In Sec. Ill we 
derive conditions for the occurrence of SDE, which are 



shown to be (almost) always met in Sec. IV. We summa- 
rize our conclusions and some possible avenues for future 
research in Sec. 0. 



II. DYNAMICS OF ONE AND TWO QUBITS 
UNDER GENERAL ERROR AGENTS 



In this section, we derive an exact analytical expression 
for the dynamics of a single qubit coupled to a Marko- 
vian environment at zero temperature for a comprehen- 
sive class of interaction Hamiltonians. We also briefly 
review two known results that will be useful throughout. 
The first one describes how the dynamics of a system of 
two qubits — each of which coupled to an independent 
environment — can be constructed from the single-body 
dynamics, whereas the second one provides a recipe to 
quantify entanglement of a two-qubit system. 



A. The Lindblad equation 

We proceed by calculating the dynamics of a single 
perturbed qubit, assuming a general interaction-picture 
Hamiltonian of the form 



Hi{t) = (A • (T)B^{t) + (A* • (T)B{t), 



(1) 



where A is an arbitrary three-dimensional vector with 
complex components, B{t) = u\,{t)BUE{t)^ f = xcr^; 4- 
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yay+zaz, and UE{t) = exp {—iHEt), with He represent- 
ing the environment Hamiltonian (hereafter, we assume 
h= 1). 

Without specifying He and B(t), the Born-Markov 
master equation for the Hamiltonian given above, at zero 
absolute temperature, T = 0, can be derived by follow- 
ing the same steps outlined in Ref. and yields the 
Lindblad equation 



dpijt ) 
dt 



= -7 



{c^Cpiit) + pi{t)C^C - 2£pj{t)C^Y (2) 



with C — X ■ cr and 

1= r ds TvE {b{0)B\s)pe{0)} . (3) 

To derive Eq. (|2[), we suppose that pe{0) is stationary, 
that is, [psiO), He] = 0, also that 



TV, 



{Bit)pEiO)}^0, 



(4) 



and that B{t) anihilates the ground state |0) of the en- 
vironment: 

m io) = 0. (5) 

Moreover, we also assume that pe{0) = |0)(0| at T = 0. 



1. Linear Dependence: Flip Hamiltonians 

Since u and v are chosen to be linearly- dependent^ we 
have w = 0, in which case it is not difficult to show that 
the solution of Eq. (0) is 



r(t) = e-47*ro + (1 - c-4^*) (ro-u)u. 



(9) 



where we have defined rg = r(0), and the unity vector u 
such that u = |u|u. 

Alternatively, it is relatively straightforward to repre- 
sent the solution above in the operator-sum form. 



p/(i)-^A',(t)pz(0)A;t(t), 

i 

with Kraus operators given by 



(10) 



^iW==VP(i)l2 and K2{t) = v^l-pW (u • cr) 



where p{t) 



(1 



-i-yt 



)/2. 



(11) 



Noticeably, the Kraus operators above are identical to 
those of the bit flip, bit-phase flip and phase flip chan- 
nels if u is chosen equal to x, y or z, respectively [ p^ , 
which motivates the name flip Hamiltonian. In fact, any 
interaction Hamiltonian arising from a linear dependent 
choice of u and v yields a type of "flip channel". The 
effect of an arbitrary flip channel on a single qubit is 
detailed in the Appendix In that appendix, we also 
illustrate the action of more general dissipative Hamilto- 
nians on the Bloch ball. 



B. Solution of the Lindblad equation 

In order to solve Eq. (||), it is convenient to use the 
decomposition A = u+iv with u, v e and |up-|-|vp = 
1 in the interaction Hamiltonian of Eq. (|l|) , in such a way 
that the Lindbladian operator can be written as 

£ = (u + iv) • cr and C'' C ^ I2 - 2w ■ a , (6) 

where w := u X V. 

Substituting the above into Eq. (|2|), and expressing 
pi{t) in terms of its Bloch vector r{t), 



I2 1.x 



(7) 



we obtain the following master equation for r{t), equiva- 
lent to Eq. (1) for pi{t), 



dv{t) 
dt 



47 I u[u-r(t)] +v[v-r(t)] +2w-r(i)| 



(8) 



In what follows, the solution of Eq. (|^) is given in two 
separate cases, depending on whether u and v arc chosen 
to be linearly-dependent or independent vectors. 



2. Linear Independence: Dissipative Hamiltonians 

In this case, since w 7^ 0, the Bloch vector can be 
expressed as 



r(t) = /(t)u + g(t)v + /i(t)w. 



(12) 



and Eq. (||) is solved with the following choices of /(t), 
g{t), and h{t): 



hit) 



and 



9 it) 



(v X w) • ro 
(w X u) • ro 



= M 

with the 2x2 matrix A4 given by 



-27* 



M = cosh(27gt) I2 + sinh(27gt) 



1 



(13) 



(14) 



, (15) 



where q = y + (u • v) and x = (|up — |vp) /2. 

From the equations above, a set of Kraus operators 
{-ft',;} for the dissipative dynamics of a single qubit can be 
obtained following a standard procedure [n9, [20(1 outlined 
here in Appendix In summary, Eqs. (pF([f5D describe 
the dynamics of a qubit coupled to a Markovian environ- 
ment at zero temperature for a broad class of interaction 
Hamiltonians parametrized by A. 
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C. Asymptotic limit 

In order to enable the study of t he e ntanglement dy- 
namics of a two-qubit system in Sec. Ill, it will be helpful 
to determine the asymptotic form of the Bloch vector of 
each individual qubit state, evolving under flip and dis- 
sipative Hamiltonians. 

In the case of flip Hamiltonians, each individual Bloch 
vector converges to Voc = (i"o • u) u [cf. Eq. (^], reveal- 
ing that the final state depends both on the interaction 
Hamiltonian (through u) and on the initial condition 
(through ro). 

In contrast, for dissipative Hamiltonians, the final state 
depends only on the interaction Hamiltonian. In fact, a 
quick inspection of Eqs. (^2|)-(|l^ shows that each indi- 
vidual Bloch vector converges to too = 2w. The inde- 
pendence of Too on the initial state conveys important 
information about the SDE phenomenon in the case of 
dissipative dynamics. 



D. From a single qubit to two independent qubits 

Since we are considering that the two qubits are inde- 
pendent and interact only locally with their reservoirs, 
their dynamics can be obtained by the set of Kraus op- 
erators that act over each single body Accordingly, 
in the case of a two-qubit initial state p{0), we have the 
following Kraus decomposition of p{t): 



p{t)= K\"{t)®KY'(t) p{0) Kl''{t)®KY'{t) 

'(16) 

where repeated indexes are to be summed over, and 
the corresponding matrices -fsTj*-"' are Kraus operators de- 
scribing each single qubit dynamics. 



,(2) 



'(1), 



E. Quantifying Entanglement 

To examine the entanglement dynamics of the den- 
sity matrixp(t) defined above, we will utilize the concur- 
rence |3|, hereafter denoted by C[p{t)]. According to 
this measure, p{t) is separable if and only if C[/ci(t)] = 0. 
On the other hand, the more entangled p{t) is, the larger 
is C[p{t)]. 

The concurrence is defined as 



CW)] 



:[0,A(i)], 



(17) 



where K{t) = h ~ I2 ~ h — h, and h > I2 > h > I4 
are the square roots of the eigenvalues of the matrix 
p{t)a^'^ p* {t)af^ , where the complex conjugate of p{t), 
denoted p*{t), is taken with respect to the computational 
basis. 

Given an initially entangled state, we note that a 
sufficient condition for the occurrence of SDE is that 
A{t 00) — henceforth denoted A^o — is a strictly 



negative number. In this case, we can guarantee that 
A(t) = for some finite time r, or in other words, that 
all the initial entanglement disappeared at r. Moreover, 
since the concurrence is known to decay monotonically 
under the Markovian approximation , in our analysis 
there will be no room for entanglement revival after r. 

In the next section, the above sufficiency criterion will 
be central in deriving more immediate conditions for the 
occurrence of SDE. Our goal is to obtain criteria that do 
not require the knowledge of the asymptotic state, but 
instead that predict the incidence of SDE from the simple 
observation of the initial states and/or the interaction 
Hamiltonian. 



III. CONDITIONS FOR SDE 

In this section, we prove two theorems that character- 
ize the occurrence of SDE in terms of quantities which 
are generally assumed to be known a priori: the initial 
state of the system and the Hamiltonian mediating the 
interaction between the system and the environment. 

Our first theorem deals with the case of flip Hamilto- 
nians, and can be seen as a generalization of a previous 
result by Huang and Zhou The second theorem, 

which treats the case of more general dissipative Hamil- 
tonians, had not yet been identified to the best of our 
knowledge. In the Appendix this second theorem is 
illustrated through the explicit analysis of some particu- 
lar cases of interest. 

Theorem 1. A necessary and sufficient condition for the 
occurrence of SDE in a system of two entangled qubits 
evolving independently under flip Hamiltonians, is that 
the matrix 

po ^ ([/(I) ® C/(2))tpo([/(i) ^ [/(2)) (18) 

has no zeroes in the diagonal. In the above, pa is the 
initial two-qubit density matrix, and J/^"-* — I2 if u'"^ = 
z, or else 



cos ip 



(")] 

^—^ sint/s^"^ 



smip 



(n) 



cos if 



(n) 



where (p^"'^ 



(19) 



(n) , r (n) (n) (n)-\ ^(n) 

iu^ ' and [u\ , , U3 J = u'- ^ 



Proof. We start by showing that, for purposes of entan- 
glement quantification, an arbitrary local flip dynamics 
of pq is equivalent to local dephasing of po- 

First, it is easy to check that [/'^"' is a unitary matrix 
such that 



(20) 



Hence, the Kraus operators of an arbitrary flip channel 
A'j"^ and A:^"^ [cf. Eq. iM)] can be related to the Kraus 
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operators of a phase-flip channel if^"-* and i^2"^ as fol- 
lows: 

^(") = [/(«)^(«)[^(n)t ^(n) ^ f;(n)_^(")[/(n)t _ 

(21) 

The two-qubit density matrix dynamics is obtained by 
plugging the initial state po and the Kraus operators 
above into Eq. (16), which becomes 



with 



Pit) 



m = 



(1) 



K 



(2) 



m 



t 



(1) 



(23) 



/or an arbitrary quhit density matrix g2, and where w'"^ 
is implicitly and exclusively defined by the Kraus opera- 
tors [i.e., Q2 is independent o/w^"^/, then 



(25) 



(22) for every 4x4 density matrix g^. 



Proof. First, notice that the linearity of Eq. (|2^) and a 
few appropriate choices of enables one to determine 
the action of the quantum channel on the operator basis 
{|M)('^|}A<,>^e{o,i}. namely 



where the Einstein convention sum is adopted. 

From the equations above, it is clear that the task of 
quantifying the entanglement of p{t) is equivalent to the 
task of quantifying the entanglement of p(<), since p{t) 
and 'p{t) are related by a simple local unitary transforma- 
tion, and thus have the same amount of entanglement. In 
other words, the problem of quantifying the entanglement 
of Po evolving under arbitrary flip channels is equivalent 
to the problem of quantifying the entanglement of the 
"rotated density matrix" evolving under dephasing 
channels. 

With this equivalence in mind, the theorem is promptly 
proved by invoking an independent result due to Huang 
and Zhu |p5| , which states that for a two-qubit system 
whose parties evolve independently under the influence of 
dephasing channels, the absence of zeros in the diagonal 
of the initial density matrix is a necessary and sufficient 
condition for the occurrence of SDE. In what follows, we 
provide an independent proof of th e sufficiency clause by 
showing that Aqo < (cf. Sec. |II E] ) , with saturation only 
if Po has at least one zero in the diagonal. 

Let r^j, with 1 < i,j < 4, be the matrix elements 
of Po- A direct computation shows that p{t — > oo) is 
simply the matrix obtained from po by zeroing the off- 
diagonal elements and keeping the diagonal ones. From 
the resulting diagonal matrix, it is easy to find that 
Aqo = — 2 min(ri_ir4_4 , r2,2''^3,3) < 0, where the inequal- 
ity follows from the fact that a positive semidefinite ma- 
trix (such as Po) can only have non-negative diagonal 
entries. Moreover, it is clear that saturation can only 
occur if ri.ir4.4 = or r2,2''3,3 = 0. □ 

In what follows we provide analogous conditions for 
the occurrence of SDE in the case of dissipative Hamilto- 
nians. Before stating and proving our conditions, it will 
be helpful to present the following lemma 

Lemma 1. If a set of Kraus operators {k\^^} induces a 
quantum channel such that 



E 



1_2 
2 



(24) 



(«)t 



12 

2 



(26) 



Now, expand g/^ in the corresponding tensor product ba- 
sis 



(27) 



where the expansion coefficients must satisfy 
J2a KPa,a,K,K = 1 in order to guarantee the unity- 
trace property of ^4. The lemma is trivially proved by 
pluggingthe expansion above into the left-hand side 
of Eq. (|2^), and by evaluating the resulting expression 
with the aid of Eq. (||). □ 

Theorem 2. A sufficient condition for the occurrence of 
SDE in a system of two entangled qubits evolving inde- 
pendently under dissipative Hamiltonians, is that the pair 
of vectors u*^"^ and v*-") specifying each Hamiltonian sat- 
isfies 



in) 



X V 



(")| 



7^ 



(28) 



for n = 1,2. Or, equivalently, none of the qubits under- 
goes a standard amplitude- damping channel. 



Proof. Recall from Sec. [I C that, in the asymptotic 
regime, a single qubit evolving under a dissipative Hamil- 
tonian reaches a state which depends only on the interac- 
tion Hamiltonian via the vector w. For two qubits, each 
of which evolving under independent dissipative Hamil- 
tonians, the asymptotic global state can be promptly ob- 
tained by invoking Eq. ( |l^ ) and the result of Lemma ^ 



I2 
2 



12 

2 



.(2) 



(29) 



which implies that any initial entanglement of po has dis- 
appeared. In what follows, we focus on determining some 
cases where such entanglement disappearance occurred 
abruptly (SDE), as opposed to asymptotically. 

Substituting the trivial identities w^"^ = |w(")|w'^"^ 

and w(") • <T = T^(")cr^M^(")^ into Eq. (|2|) [in which the 
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unitary matrix W^") is defined in analogy to Eqs. ( |l9| ) 
and (I20I)], one finds that 



witli 



I2 
2 



|w(i)|a. 



I2 
2 



Iw(2)|a, 



(30) 



(31) 



Tlrus, without loss of generality, the task of quantifying 
the entanglement of p^o can be replaced with the task 
of quantifying the entanglement of the diagonal density 
matrix poo, for which one can easily compute 

Aoo = -^V^[l-|2w(i)|2] [l-|2w(2)|2] <0. (32) 

In the above, the inequality is trivially satisfied with sat- 
uration if jw*^")! = 1 /2 for 71 = 1 or n = 2. As noted 
before (cf. Sec. p^I E| ), a strictly negative value of Aqo 
indicates the disappearance of entanglement in a finite 
time interval, leading to the sufficient condition for SDE 
given by Eq. (^). Remarkably, this is not a necessary 
condition since Aqq may have vanished at a finite time 
and remained zero ever since (cf . Appendix |c|) . 

Finally, we note that if |w(")| = 1/2, then the n-th 
qubit collapses to a pure state, regardless of the initial 
preparation. As seen in Appendix this is precisely 
what is accomplished by a standard amplitude damping 
channel. Therefore, with this only exception, SDE occurs 
for all possible choices of dissipative Hamiltonians. □ 



We now show that the identified dynamical models for 
which SDE does not necessarily occur form sets of di- 
mension 3. Let 



SIf = {(u,v) : V = au,a € M}, (35) 
r!AZ5 =f^n{(u,v) : |uxv| = i}, (36) 



where ^Ip is equivalent to the set of flip Hamiltonians, 
and ^AD is equivalent to the set of Hamiltonians pro- 
ducing standard amplitude damping channels (cf. Ap- 
pendix A 2). 

A straightforward computation shows that any ele- 
ment oi Vtp can be parametrized as in Eq. ( ^^ with 
i? = (1 + a^)"^/^, 0' = 6 and <f>' = cf), in such a way 
that two free parameters are eliminated. As a result, we 
have dimQp = 3. 

Similarly, one can easily show that the elements of ^ad 
can be parametrized as in Eq. (11) with R ^ 1/^2, and 
one of the angles fixed by the solution of the transcen- 
dental equation u • v = 0. As a result, we also have 
dim flAD = 3. 

Since the volume of f2i? and ^Iad with respect to the 
standard Lebesgue measure in vanishes, the theorem 
is proved. Pictorially, one should consider the analo- 
gous lower dimensional scenario of a curve lying inside 
a solid block: it is impossible to cover the whole block 
with a countable number of curves, so the curve has zero- 
measure with respect to the block. The same argument 
holds for higher dimensions. 

□ 



IV. UNIVERSALITY OF SDE 



With the aid of the conditions for SDE derived in the 
preceding section, we are now ready to prove our main 
result. 

Theorem 3. The set of dynamics for which SDE is not 
guaranteed to occur has zero-measure in the set of all 
possible dynamics. In other words, SDE is universal. 

Proof. For a single qubit, the set of all possible dynamics 
generated by the Hamiltonian of Eq. (|l]) is equivalent to 
the set 

f} = {(u,v) : u,veRMu|2-f |v|2 = 1}, (33) 

thanks to the parametrization A = u + zv with |A| = 1. 

The dimension of a set is defined as the number of 
independent parameters needed to describe a point in 
the set, in which case dimfi = 5, since any element of 

can be written in terms of the parameters R E [0, 1], 
9, 0' e [0, 2tt] and 0, 0' G [0, tt] as follows: 

u = i? (cos 6 cos X + sin 6 cos (f> y + sin z) (34a) 

V = -v/l — i?2 (cos 9' cos (/)' 'X. + sin 9' cos (f>' y + sin (/>' z) 

(34b) 



V. CONCLUDING REMARKS 

As far as nature does not present any preferential 
way non-interacting qubits can couple to their assumed- 
Mar kovian environments, SDE is almost-surely universal 
even at zero temperature. The universality achieved in 
the present work is conditional to the stated assumptions 
which, by their experimental banality, are to be com- 
monsensically deemed ubiquitous, except for the zero- 
temperature restriction. Nonetheless, recent theoretical 
investigations demonstrate that, at finite temperatures, a 
relevant set of two-qubit entangled states, robust against 
SDE at zero temperature, get disentangled at a finite 
time ||2^. The results of those investigations, although 
also obtained under the assumption of non-interacting 
qubits coupled to Mar kovian environments, give credence 
to our intuitive notion that if SDE occurs at low temper- 
atures, it necessarily occurs at higher temperatures. Of 
course, a thorough demonstration of universality is be- 
yond the scope of the present work and approaching the 
case of non-Markovian environments at finite tempera- 
ture constitutes a natural extension of current research. 
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APPENDIX A: IN-BUILT STANDARD 
DYNAMICS 

In this section, wc illustrate the generality of our noise 
model by showing that it encompasses various standard 
single qubit dynamics. This section in divided in two 
parts: In the first, we show that flip Hamiltonians lead 
to "flip channels" such as bit flip, phase flip and bit- 
phase flip. In the second part, we show that dissipative 
Hamiltonians produce dynamics represented by channels 
such as depolarizing and amplitude-damping. 



1. Flip Channels 

In Eq. (H), we found that the Bloch vector of a single 
qubit evolving under a flip Hamiltonian always points 
towards a direction u S R'^ specified by the Hamiltonian. 
Now let U"'" and u'~ be orthogonal directions to u such 
that {u, u-*", u'"} forms an orthonormal basis for M'^. We 
find that the component-wise evolution of r{t) satisfies 



As it should be expect ed, th is observation is consistent 



r{t) • u = ro • u 
r{t) ■ u-L = e-^''' (ro • u^) 



(ro • ii^) 



(Ala) 
(Alb) 
(Ale) 



which shows that while the u-component remains un- 
changed, the orthogonal components exponentially de- 
cay in time at a rate A-f. The effect of such a channel 
on the whole Bloch sphere can be visualized in Fig. |l| 
for two time-steps and an arbitrarily chosen direction u, 
designated by the arrow. 




FIG. 1: (Color online) In the first Bloch sphere, yt — 0.3, 
while in the second, — 0.7. The arrow designates the di- 
rection u. 

From the figure, it is clear that, asymptotically, the 
Bloch ball collapses to a diameter of the original sphere. 



with the result of Sec. II C, which establishes that r, 



is not completely specified by this type of interaction 
Hamiltonian, but also by Tq. That is because different 
initial states can collapse to different points of a particu- 
lar diameter (specified by u). Furthermore, it is obvious 
from Fig. [l] that if u were chosen to be one of the Carte- 
sian versors x, y, or z, then the resulting channel would 
be, respectively, a bit flip, a bit-phase flip, or a phase flip. 

Notice that Fig. Q illustrates the most general dynam- 
ics that can arise from a flip Hamiltonian. Nevertheless, 
this is certainly not the most general dynamics that there 
is. For example, in the present case the Bloch ellipsoid 
will always be centered at the origin [unital channel), 
and will always touch the Bloch sphere in two antipodal 
points, meaning that there is always a pair of orthogo- 
nal pure states that remain pure even after coupling to 
the environment. Needless to say, these features are not 
generally observed in quantum dynamics. 

In the next section, we show that dissipative Hamilto- 
nians give rise to a whole new class of channels, including 
non-unital channels such as the amplitude-damping, and 
unital channels that globally drop purity, such as the de- 
polarizing channel. 



2. Amplitude-Damping and Depolarizing Channels 

Here, we take A = u -I- iv with u and v linearly- 
independent vectors. In particular, we choose 



u = x sm (p cos t/ -h z cos ( 
V = — y sin sin 9 



(A2) 
(A3) 



in which case the normalization condition |A| = 1 is ful- 
filled for any 9 G [0, 2tt] and cj) e [0, tt] 

Now wc look at the time evolution of the components of 
the resulting Bloch vector of Eq. (p^). Note that u and 
V are orthogonal vectors, in such a way that {u, v,w} 
forms an orthonormal basis for M"^, where wc have 



w u X V = sin 9 sin (f) (x cos — 2 cos 9 sin i 



(A4) 



and w := w/|w|. Substituting the above expressions for 
u, V and w into Eqs. (^3|)-(^5|), the dynamics of the u, 
V, and w components of r(t) is given by 



r(t) •u==e-27-t(rp.u) 



r(t) • V 



-27+t 



(ro • v) 



v{t) ■ w = e"'*'^*(ro • w) + 2(1 - e 
where we have defined 

27± := 27 ± [cos(26l) sin^ cj) + cos^ 



-47* 



)|w| 



1 1/2 



(A5a) 
(A5b) 
(A5c) 

(A6) 



By comparing Eqs. (|A5| ) and (Al), we notice two fun- 
damental differences in the resulting dynamics of the 
Bloch vector: First, in the present case, there is no longer 
an invariant component as the one implied by Eq. ( Ala| ). 
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Second, as t ime goes by, the second summand on the rhs 
of Eq. ( |A5c| ) shifts the origin of the Bloch sphere along 
the w-direction by an amount proportional to |w|, hence 
giving rise to a non-unital channel. 

In this regime, it is thus clear that unital channels can 
only be obtained if |w| = 0, or equivalently, if 9 and (j) sat- 
isfy cither (i) sin = or (ii) sin = or (iii) cos </> = 
and sinO ~ 0. However, in contrast to the previous sec- 
tion, the unital channels arising here require the lengths 
of all the Bloch vector components to exponentially de- 
cay with time, asymptotically collapsing the Bloch ball 
into a point corresponding to the maximally mixed state 
I2/2. A unital channel of this form is known as a depo- 
larizing channel, and is illustrated in Fig. || for two time 
steps. 




FIG. 2: (Color online) Depolarizing channel. 9 = tf) = 0. In 
the first Bloch sphere, jt = 0.3, while in the second, — 0.7 

As noted before, for |w| 7^ a family of non-unital 
channels arises. A well-known representative of this fam- 
ily is the standard-amplitude-damping channel, occur- 
ring whenever \9\ = tt/A and (f> = tt/2. The action of 
this channel (with 9 = +7r/4) on the Bloch sphere is 
illustrated in Fig. || for two time steps. 




FIG. 3: (Color online) Amplitude-damping. 9 = n/A and 
= 7r/2. In the first Bloch sphere, 'yt — 0.3, while in the 
second, "/t = 0.7 

In analogy to the case of depolarizing channels. Fig. ^ 
illustrates that the Bloch ball asymptotically collapses 
to a single point. The difference between the two cases, 
though, is that for amplitude-damping this point corre- 
sponds to a pure state, as opposed to the maximally- 
mixed state in the case of the depolarizing channel. 
Clearly, the property of mapping the Bloch ball into a 



single point is nothing but a direct consequence of the 
limit establishing that Too is independent of Tq for this 
type of interaction Hamiltonian (cf. Sec. |II C| ). 

Finally, it should be noted that generalized amplitude- 
damping channels can also be obtained by varying the 
parameters 9 and (j). In those more general cases, the 
Bloch sphere can collapse to states that are neither pure 
nor maximally mixed, as revealed by a quick inspection 
of the magnitude of the vector w of Eq. (A4). 



APPENDIX B: KRAUS OPERATORS 

In this appendix, we show how to construct any set 
of Kraus operators for a channel acting on qudit states. 
Besides being of independent interest, this result is use- 
ful to derive a set of Kraus operators to represent the 
dynamics of a single qubit evolving under a dissipative 
Hamiltonian (cf. Sec. II B 2 ). 

Before proving the main result, let us recall a few useful 
facts from the theory of quantum channels p7| . 
Let Cl*^* denote the set of all quantum channels C acting 
on d-dimensional density matrices, and let denote 
the set of all d^-dimensional positive semidefinite matri- 
ces £ such that Tr2 € — Id- It was independently shown 
in Refs. that C§^* is isomorphic to C^l*, with an 

isomorphism law — sometimes referred to as the Choi- 
Jamiolkowski isomorphism — given by 



E \j){k\(SC{\j){k\) 



(Bl) 



where is an orthonormal basis for the Hilbert 



space of dimension d. 
verified M, 



Conversely, as it can be easily 



C(p) -Tn [{p'®td)€] 



(B2) 



where the transposition is taken with respect to the basis 

The crucial point to bear in mind is that the knowl- 
edge of a quantum channel C is entirely equivalent to the 
knowledge of the corresponding matrix £, henceforth de- 
nominated the Choi matrix of C |31 . In what follows. 



we characterize the Kraus operators of C in terms of a 
decomposition of C 

Lemma 2 (Verstraete Salgado Let p{t) = 

Ct{po) be the density matrix of a qudit at time t, where 
Ct is a quantum channel and po = p(0). Moreover, 
let £(f) be the Choi matrix of Ct, and consider any d^- 
dimensional matrix S{t) such that S{t)S{ty = 2^(0- ^ 
Kraus decomposition for p[t) is 

p{t) = E '^^(^)] ^'(^)]^ 

i 

where Si{t) is the i-th column of S{t), and mat Si{t) is the 
dx d matrix whose n-th column is formed by the elements 
of Si with indexes ranging from {n — l)d + 1 to nd. 
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Remark 1. The operation mat is simply the inverse of 
the operation vec i^^, which vectorizes a matrix by stack- 
ing its columns. For example: 



vec 



a c 
b d 



id mat 



b d 



(B4) 

Remark 2. The Choi matrix of the quantum opera- 
tions described in Sec. II B I can be easily computed via 
Eq. (^) by noticing t/iaTCljO) (0|j_an(i C(|l) (1|) can be 



directly obtained from Eqs. (|12[)-(|15|) by making Tq — z 
and Tq = —z, respectively. The matrices C(|0)(1|) and 
C(|1)(0|), in turn, can be obtained by exploiting the lin- 
earity of C , explicitly 



2C(|0)(1|) = 
C(|+)(+|)-C(|-)(-|) 



iC{\+i){+i\)-iC{\-i){-i\) 
(B5) 



-i\)+iC{\-t){-z\) 

(B6) 



2C(|1)(0|) = 

c(l+)(+l)-c(|-)(-|)-*c(|+*) 



where C{\±){±\) and C{\±i){±i\) follow from Eqs. (|2|)- 
( p^ ) by making ro ~ ±x and Tq ~ ±y, respectively. 

Remark 3. Since €{t) is a positive semidefinite matrix, 
it is always possible to find S{t) such that S{t)S{t)'^ = 
fL{t) (e.g. Cholesky decomposition). In fact, there are 
(infinitely) many such choices of S{t), which implies in 
(infinitely) many Kraus decompositions of p{t) 



Proof. The proof is a straightforward apphcation of the 
foUowing easy-to-check matrix identity: 

ABC^ = Tn [{B'^ (E) Id) (vec A) (vec C)^] (B7) 

for A,B,C arbitrary d-dimcnsional matrices. In partic- 
ular, take A = C = mat Si and B ~ po, then 



[mat 5,] po [mat 5*,]^ = Tn {pi ® 1^) S^Sj 



(B8) 



where we have used that v ec [m at Sj] = 5*^. Summing 
over i on both sides of Eq. (B8) and exploiting the lin- 
earity of the partial trace, the desired result follows: 



[mat Si] Po [mat Si]'' ~ Tri 



[Po 



= Tn [{pI ® Id) ^^t] 
= Tn [{pI ® Id) £] 
= Ct{po) 

= Pit) (B9) 



□ 



APPENDIX C: WATCHING SDE 

As an illustration of the result of Theorem ||, in this 
appendix we aim to visualize the occurrence of the SDE 
for some instances of two-qubit entangled states, whose 
parties evolve independently under identical dissipative 
Hamiltonians. 

For the initial conditions, we consider the states 



a-lTi>+/3-liT>, 
a+lTT>+/3+Ui), 



(Cl) 
(C2) 



where, for simplicity, we suppose a±, (3± £ M. 

For the Hamiltonians, instead of the usual absorption 
coupling agents a± = {ux ± iay)/^/2 — or, in our nota- 
tion, u*^"^ = x/-\/2 and vf^"^ = y/v^ — we consider the 
more general dissipative couplings given by 



and 



ysmi 



(C3) 



for e G [0,27r] and n = 1,2. 

With the initial states and Hamiltonians thus specified, 
the dynamical equation (|l^), and the Kraus operators, 
it is possible to obtain the density-matrix dynamics for 
continuous time t. Furthermore, with a simple numerical 
procedure, the function K{t) can also be computed in 
each case, and is plotted in Figs. ^ and |[ 
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FIG. 4: (Color online) I^^{t) for the initial conditions given 
by ja+j^ = 0.2 for the red (solid) line, |q:+|^ = 0.5 for the 
black (double-dotted) line, and |q:+|^ = 0.8 for the blue (dot- 
ted) line. In (a) and (b) we use the interaction Hamiltonian 
as defined by Eq. (C2), with 6 = — 7r/4 and 6 — —it/5, re- 
spectively. 



Assuming 7 
given by Eq. 
= 0.2, la 



1 



in Eq. (y) and the initial condition 

( |C2| ), we illustrate, in Fig. |, A+(t)_for 



0.5, and 



0.8. In Fig. 1(a) 



we consider t he u sual amplitude-damping, supposing 6 ~ 
— 7r/4 in Eq. (|C^), and we observe that the entanglement 
suddenly vanishes for \a+ [^ = 0.8. In fact, as it is already 
known in these circumstances ^, |l^, ij, |l5[ p^ , 
SDE occurs whenever \a_+\'^ > 0.5. In Fig. E(6) we con- 
sider 6 = — 7r/5 in Eq. (C3). In this case there is no SDF 
subspace for the initial conditions of Eq. (K 
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FIG. 5: (Color online) A~(t) for the initial conditions given 
by = 0.2 for the red (solid) line and \q—^ — 0.5 for the 

black (double-dotted) line (for this class of initial conditions, 
|q_P = 0.2 and |a-P — 0.8 give identical dynamics). In 
(a) and (b) we use the interaction Hamiltonian as defined by 
Eq. (C3), with 6 = — 7r/4 and 9 = —it/5, respectively. 



In Fiff. I wc consider the initial conditions given by 
Eq. (|Cl[) and keep all other variables as in Fig. |[ In 
this case, as illustrated in Fig. ^(a), we have a SDF 
subspace for the usual amplitude-damping {9 = — 7r/4). 
Otherwise, as Fig. ||(6) shows, SDE is present for all 
initial conditions when 8 = — 7r/5. 



Finally, we note that the observations above are fully 
consistent with the statement of Theorem ||. The 
fact that SDE can occur even in the case of standard 
amplitude-damping channel proves that the condition of 
the theorem is not a necessary one. Moreover, the occur- 
rence of SDE whenever a channel slightly different from 
the standard amplitude-damping is applied, supports the 
sufficiency clause. 
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